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The relationship between the volume and surface energy coefiicients in the liquid drop A~^^^ 
expansion of nuclear masses is discussed. The volume and surface coefiicients share the same physical 
origin and their physical connection is used to extend the expansion with a curvature term. A 
possible generalization of the Wigner term is also suggested. This connection between coefficients 
is used to fit the experimental nuclear masses. The excellent fit obtained with a smaller number of 
parameters validates the assumed physical connections and the usefulness of the curvature term. 



Introduction. Nuclear masses and their dependence 
on atomic and mass number gave essential information 
about the nature of nuclear forces. They also led to the 
formulation of the liquid drop model, arguably the most 
precise and easily interpretable description of the masses 
themselves 

Empirical trends and scientific intuition led to the for- 
mulation of the liquid drop model. In its traditional form, 
the liquid drop model approximates the binding energy 
of a given nucleus of mass number A and charge Z as [l| : 



(A-2Z)2 S 



The five terms in this equation are associated with five 
independent aspects of nuclei expected to affect the bind- 
ing energy. These aspects are the nuclear volume, sur- 
face. Coulomb repulsion, proton-neutron asymmetry, and 
pairing. A fit of this equation to nuclear masses gives 
the coefficients and reproduces the experimental values 
to within 1% or ^10 MeV for heavy nuclei. 

This is an outstanding result that attests to the pro- 
found physical content of the overall equation and to the 
interpretation of its individual terms. The residual 1% 
discrepancy is due to shell structure. The shell correc- 
tions, evaluated according to the Strutinsky procedure 
01 and grafted onto the liquid drop model, permit an ac- 
curate evaluation of nuclear masses and fission barriers 
to within 1-2 MeV This hybrid approach remains 

to this day the yet unmatched paragon for more sophis- 
ticated models such as Hartree-Fock-Bogoliubov :7i, 8|. 

Motivated by this early success of the liquid drop 
model, many additional terms have been suggested, each 
with its own physical interpretation. An example of this 



is found in Myers and Swiatecki [j^ who suggested: 
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where I~A — 2Z. The main difference between Eq. ([T]) 
and Eq. ([2]) is the extension of the neutron-proton asym- 
metry to the surface energy term. Also, a term linear in 
|/| was introduced. 

The further physical insight found in writing the asym- 
metry energy as in Eq. ([2]) is the connection it implies 
between volume and surface energies. The authors ar- 
gued that the change of the volume energy due to the 
neutron-proton asymmetry I should be refiected in the 
surface energy of the system as well, though stating that 
this was done without empirical evidence [3]. The nat- 
ural implication is that the surface and volume energies 
are related through their common origin. 

A term linear in |/| was originally suggested by Wigner 
in considering the exchange force of nucleons [9]. An 
empirical observation of such a dependency in the masses 
was reported by Myers and Swiatecki, hence its addition 
in the above equation [3|. 

The last alteration was the addition of a term propor- 
tional to Z'^ /A, which is to account for the difference in 
the Coulomb energy due to the diffuse nuclear surface. 
Again, this trend was not observed in the nuclear data, 
but was added because it would be anticipated. We will 
not discuss the role of this term in the nuclear binding 
since it requires an additional fitting parameter to de- 
scribe what is ambiguous in the data. 

In this paper we extend and generalize the insights dis- 
cussed above. In particular, we argue that the relation 
between volume energy and the surface energy is strong 
enough that their coefficients should not be taken as in- 
dependent variables. Furthermore, we discuss the need of 
a third term arising from the same physics, proportional 
to A'^/'^, in order to create a consistent physical picture 
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of the nuclear binding energy. Finally, a linear term in 
|/| is naturally introduced when treating the asymmetry 
term as the expectation value of the isospin, T^. 

A revised liquid drop model is fit to the experimental 
binding energies of the nuclides to test these considera- 
tions. The importance of the revised terms is assessed by 
comparing the fit of the original and revised models. 

The Liquid Drop Formula: A Truncated Series 
Expansion. The first term of Eq. ^ is aptly called the 
volume term. Its proportionality to A indicates saturat- 
ing forces leading to constant density and binding energy 
per nucleon. The obvious similarity to molecular fluids 
led naturally to the introduction of the second term, the 
surface term. Its proportionality to A^^^ speaks to the 
lack of saturation on the nuclear surface, whose area, 
through the constant density of the fluid, should indeed 
be proportional to A^/^. Progressing along the same line, 
it was widely appreciated that the surface term is a finite 
size correction and that additional terms in the expansion 
might be needed, such as a curvature term. 

Generally, we can think of a generalized liquid drop 
formula as a rapidly converging series expansion in pow- 
ers of A~'^/'^ , known as the leptodermous expansion [lo| : 

EB^-a^A + asA^/^ +arA^/'^ + ... (3) 

It is left to be determined how many terms in the expan- 
sion are necessary to describe the physics of the nuclear 
system. The incorporation of a curvature term, with its 
coefficient ar, proportional to A'^/^ is almost demanded 
by the truly small size of nuclei {A < 300) compared to 
the size of the drops typically considered in molecular flu- 
ids, such as aerosols, where A > 10^. Higher order terms 
also may be of importance due to the small size of nu- 
clei, but would be intractable without an understanding 
of the lower order curvature term. 

The role of the curvature term in nuclear systems was 
considered only recently and has yielded ambiguous re- 
sults 3- 6]. The increased number of parameters and the 
ability of the traditional liquid drop formula without cur- 
vature to fit the data made the problem of identifying the 
magnitude of this term rather difficult. We believe that 
it is possible to shed additional light on this subject by 
considering the physical origin of the various terms. 

Volume and surface terms both arise from the same 
physical property of nuclear forces: saturation, and the 
lack thereof. Thus, surface and volume terms should be 
related to one another, being themselves different effects 
of the same cause. Furthermore, the experimental sur- 
face and volume coefficients turn out to be approximately 
equal. Is this an accident or could they possibly be equal? 

To answer this question, consider a system of small 
sticky cubes used to build larger, composite cubes. These 
cubes interact only when in direct contact. The system is 
characterized by some bond strength, e, when two faces 
are touching. The energy of a cube of A constituents is 
equal to a volume energy minus a surface energy, just as 
in the nuclear case. Counting the number of bonds in a 



TABLE I. Fits of the nuclear masses with Eq. (fTTj) using dif- 
ferent mass ranges and setting ar—O. All the parameters in 
units of MeV. The value in the parentheses is the uncertainty 
in the last digit. 



Masses 


a„ 


as 


k 




5 


50-100 
100-150 
150-200 
200-250 


15.39(4) 
15.39(2) 
15.11(2) 
15.18(6) 


16.81(10) 
16.68(7) 
15.66(8) 
15.7(2) 


1.742(7) 
1.771(3) 
1.748(3) 
1.768(5) 


0.686(3) 
0.6917(14) 
0.6760(12) 

0.686(3) 


10.3(5) 
12.4(3) 
13.5(3) 
13.3(4) 



cube of size A reveals: 

Thus, in this model the volume and surface energy coef- 
ficients are exactly equal with ay^as=3e. Even though 
this is a simplified model in comparison to a nucleus, it 
exemplifies the fact that the volume and surface terms 
are strongly connected. This insight motivates setting 
ay=as without any loss of information. 

One difference between this simple model system and 
a nucleus is the diffuseness of the nuclear surface. What 
effect does a diffuse surface have on the binding energy 
of a drop? Since the volume energy is a property of the 
bulk system it would remain unchanged. The fact that 
the system naturally becomes diffuse means that it gains 
a larger binding energy in doing so. The surface energy 
would then be lowered in comparison to the sharp surface 
system. This implies that the surface energy coefficient 
should be equal to or smaller than the volume energy 
coefficient, contrary to what is observed in traditional 
liquid drop fits to the nuclear masses. 

As the system is made smaller, more terms in the lep- 
todermous expansion may be needed to properly predict 
binding energies. If one were to fit the expansion with an 
insufficient number of terms, what ailments would be ob- 
served? The terms included in the equation would have 
to change from their nominal values to accommodate the 
lack of higher order terms. Also, the deviation from the 
nominal value would be worse for smaller masses, where 
the higher order terms are more important. 

As an example, consider nuclear binding energies in 
various mass ranges. Each mass range can be fit with 
Eq. (fTTj) . using the fitting procedure that is described in 
the following section. Table U shows the results of such 
an exercise. Most terms do not vary systematically as 
the mass range is changed, their variation being of the 
order of 1%. The exceptions are the surface energy and 
the pairing energy. The pairing energy is of unrelated 
physics and is not discussed here. The surface energy 
coefficient decreases as the mass range is incremented. 
This trend indicates that the A^/^ term is not sufficient 
in describing the lack of saturation in the system. As the 
masses used in the fit increase, the surface term tends to 
the value of the volume coefficient. Hence, both the need 
of a curvature term and setting ay=as are motivated. 



3 



vacuum 




FIG. 1. Schematic representation of the surface energy. The 
system on the left represents a flat surface of an infinite liquid 
and the system on the right is a finite liquid drop. The surface 
area of a constituent of radius r„ exposed on the surface of a 
liquid drop of radius R is more than that of a surface particle 
at the fiat surface, as emphasized by the bold curve. 

Now that we see the need for a curvature term, one 
may ask: What is its origin? To answer, let us consider 
a simple liquid with spherical molecules of radius r„ . As 
shown on the left side of Fig. [U on a flat liquid surface 
the molecules should protrude half way on average, los- 
ing half of their binding energy. If the liquid surface is 
curved like that of a sphere, as in the right side of Fig. [T] 
with a drop radius of R, the molecules protrude more, 
losing additional binding. Thus, the curvature and sur- 
face terms arise from the same physical effect and their 
coefficients should be related. 

In order to obtain a quantitative estimate of this effect, 
consider a model that is geometric in nature. The surface 
energy is considered to be proportional to the protrud- 
ing surface area of a constituent residing on the surface 
times the number present on the surface. As a function 
of nuclear radius R, the resulting exposed surface area S 
of a constituent on the surface is: 

with the limiting case of a planar system, S=2Trr^. The 
number of particles on the nuclear surface is proportional 
to A^^^. The overall surface energy is then: 




Since nuclei exhibit a saturation density, the nuclear ra- 
dius is approximated as R=roA^/^ , with ro being a con- 
stant. Inserting this relation into Eq. ^ yields: 

E,^asA'/^ + as^A^/\ (7) 

Here we identify the usual surface term proportional to 
A^^^. Furthermore, we notice a curvature term propor- 
tional to A^/"^ with a coefficient that is dependent on the 
surface energy coefficient and the ratio of the "molecule" 
radius to ro, which is directly related to the saturation 
density. 

The above equation is reminiscent of the Tolman cor- 
rection to the surface tension [ll|. This term can be 
interpreted as the Tolman correction for the nuclear sys- 
tem in its ground state. 



Naturally, deviations from sphericity of the molecules 
would involve a (temperature dependent) reorientation 
on the surface. This would alter the simple relationship 
between volume, surface and curvature energies. We will 
limit the discussion to the case of an isotropic force for 
the model presented here. 

We may check the model further by putting experi- 
mental values into Eq. ([7]). Taking tq ~ 1.2 fm Q and 
the radius of a free nucleon to be r„ ~ 0.9 fm yields: 

0.9 3 

These geometric arguments thus give a first order approx- 
imation as to the sign and magnitude of the curvature 
term. Other aspects might influence the actual value of 
the curvature term in the nuclear system, but it would be 
notable if a proper fit to nuclear masses were to produce 
a value close to the above estimate. 

To further appreciate the significance of the relation 
between volume, surface and curvature energies, consider 
the following. What information is gained in knowing 
the leptodermous expansion for an arbitrary liquid in its 
ground state? 

First, consider the volume energy. The volume energy 
gives no information of the internal structure of the sys- 
tem. It is just the scale which sets the size of the rest of 
the terms in the leptodermous expansion. 

Now a measurement of the system's surface energy is 
made. The particle density of the system can be deduced 
by comparing the surface and volume energies. This is 
done by anticipating that the two coefficients will be the 
same in terms of A and A^^^, respectively. Avogadro's 
number could thus be inferred. 

Finally, the curvature energy is determined and from it 
the size of a single particle in the liquid can be estimated. 
This is shown in Eq. ([?])• 

Here we see how the hierarchy of terms in the lepto- 
dermous expansion can be related to the internal struc- 
ture of a ffuid. Even though this exercise is pedagogi- 
cal in nature, it demonstrates the physical significance of 
each term. It could have allowed Democritus to prove his 
atomic theory, had he been inclined to do so. 

Nuclear Mass Fit Results. We use a set of 2076 
masses, corrected for microscopic effects according to 
Moller et al. [13]. These microscopic corrections account 
for the shell effects along with the effects associated with 
nuclear deformation. The masses considered in the fits 
correspond to nuclear masses from reference [l3| with 
N > 7, Z > 7, and with experimental uncertainties less 
than 150 keV. The lower limit of neutron and proton 
numbers is chosen to ensure that the included nuclei are 
large enough to be considered as liquid drops. The re- 
striction on the experimental uncertainties is not only 
due to the error of the mass, but also to the reliability 
of the shell correction for masses far away from stability. 
The binding energy, Eb, of each nucleus is defined as: 

Eb {A, Z) = Znip + {A-Z)mn- M {A,Z) + A.heii {A,Z), 

(9) 
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TABLE II. Fits from the four different mass equations as described in the text. All parameters are in units of MeV. The value 
in the parentheses is the uncertainty in the last digit. 



Fit 


a„ 


as 




k 




S 


r„ (fm) 


x' 


A 


15.597(7) 


17.32(2) 


ar~0 


1.8048(9) 


0.7060(4) 


11.4(2) 




0.58 


B 


14.843(3) 


Ov—as 


ar—O 


1.7196(16) 


0.6585(4) 


10.1(6) 




4.24 


C 


15.25(3) 


15.17(17) 


3.8(3) 


1.779(2) 


0.6932(11) 


11.3(2) 


0.60(5) 


0.54 


D 


15.264(4) 




3.60(3) 


1.7805(8) 


0.6938(3) 


11.3(2) 


0.566(5) 


0.54 



with rup and m„ being the mass of a proton and neutron, 
respectively, M is the experimental mass of the nucleus, 
and Asheii is the shell correction. The liquid drop formula 
is fit to this binding energy with each nucleus given an 
equal weight. The mean square deviation of the fit is 
used to evaluate its goodness: 



(ex) T7^{th)^2 



^ N 
We use the following liquid drop formula: 



Eb =(-a„A + asA^/^ + a^A^/^) ( 1 - fc 
Z{Z-1) ^ 6 



\mi\ 



(11) 



where we insert the mass asymmetry dependence I~A — 
2Z both in the volume and surface terms according to 
Myers and Swiatecki Q. If the mass asymmetry term 
is interpreted as an "isospin" dependence, the term lin- 
ear with P should be treated as T^, with T=\I\/2. 
This "isospin" presents itself as the square T^, which we 
rewrite (with a possibly unjustified quantal sensitivity) 
as {T^)^T{T + 1) = |/|(|/| + 2)/4. This introduces a lin- 
ear term in |/| without the addition of a new parameter, 
as opposed to a freely varying Wigner term [Qi]. 
The following fits are performed: 

A. and vary independently without a curvature 
term. 

B. Same as above, but forcing a^j—ag. 

C. ay and as vary independently with a curvature 
term. 

D. Same as above, but forcing ay—a^. 

The Coulomb, mass asymmetry and pairing coefficients 
are left as free parameters in all of the above fits. The 
results are shown in Table and are discussed below. 
Fig. [2] shows plots of the residual masses of the fits, the 
exact binding energy with shell corrections included mi- 
nus the binding energy predicted from the fitted formula. 

Comparing fits A and B shows that setting a„=as with- 
out the curvature term does not ameliorate the situation. 
Quite to the contrary, the value is 8 times larger and 
the plot of the residual masses shows clear deviations. 
Left without constraint, the surface term incorporates 
the curvature effects and becomes larger. 




ay =15.60; as =17.32; Or =0.00 
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FIG. 2. The residual mass from the corresponding fits. The 
label in the top left corner of each plot corresponds to the 
fits listed in the text. The connected lines represent chains of 
isotopes. 



Comparing fits A and C, we observe that the introduc- 
tion of the curvature term as a free parameter improves 
the resulting fits as expected. But is the value of 
physically meaningful and how does it compare to the 
expectations of the geometric model? 

Rearranging Eq. ^ gives the radius of the nucleon as 



2.4 



ar 



(12) 



in units of fm. Using this equation, the nucleon radius 
is found to be 0.60(5 ) fm , smaller than the experimen- 
tal value of 0.84 fm [12|- This size of deviation is not 
unexpected from the crude approximations used, and it 
remains impressive that both the sign and relative mag- 
nitude are predicted. Furthermore, the surface energy 
coefficient moves within error of the volume energy coef- 
ficient. The other parameters change within 2% between 
the two fits, showing consistent results. 

By forcing a^—as with the presence of the curvature 
correction, as in fit D, the changes by a fraction of 
a percent. Also, the parameters not associated with the 
saturating nuclear force are left unchanged. Thus, no 
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TABLE III. Fits of the nuclear masses to the liquid drop 
model using different isospin dependencies. The first sets 
(/^) = |/|(|/| + 2), where as the second represents a fit to 
(/2) = |/j(|/| + x). All parameters are in units of MeV. The 
value in the parentheses is the uncertainty in the last digit. 

flt, Or k X ac 5 

15.264(4) 3.60(3) 1.7805(8) 2 0.6938(3) ff.3(2) 0.54 
15.247(4) 3.76(3) 1.7944(10) 1.51(3) 0.69f3(3) ff.3(2) 0.46 



physics is lost with setting a^^ag. 

Without taking into account the curvature term, the 
volume and surface parameters will tend to be irreconcil- 
ably different to be considered equal. This explains the 
reason why the two terms have previously been treated as 
independent values. The addition of the curvature term 
corrects this discrepancy, and it is found that the surface 
and volume energies are close to being equal, giving no 
visible difference in the fitting of the experimental data. 

Another fit was performed using (/■^) = |/|(|/|-|-a:), with 
the added fit parameter x. This addition is equivalent to 
introducing an adjustable Wigner term linear in isospin. 
Table IIIII shows the fit with and without letting x vary. 
None of the other fit parameters change substantially. As 
for X itself, it is found to be 1.51(3), which slightly lowers 
the of the fit. When written in the form presented by 
Myers and Swiatecki 3] , this corresponds to a congruence 
energy of 41.3(8) MeV, which agrees with the value 42 
MeV which they report. With most of the parameters 
changing less than 1%, the same physics is still captured 
by setting (/2) = |/|(|/| + 2). 

Implications of the curvature term. The existence 
of a curvature energy, especially important in light nuclei, 
may imply effects hitherto undiscovered. We give here 
two examples. 

The curvature of the surface in the nuclear deformation 
landscape, and in particular at the fission saddle point, 
exhibits large variations going from positive to negative. 
Therefore, the prediction of fission saddle point configu- 
rations and masses will be affected by the presence of a 
curvature term, which will acquire a tensorial form. 

The fragment distribution predicted by the Fisher 
model |15|] is dependent on the surface energy of the clus- 
ters. The theory uses a term proportional to A" for this 
purpose. Since the fragment yields are weighted heavily 
towards lighter fragments away from the critical temper- 
ature, the introduction of a curvature term would seem 
imperative. Thus, the curvature term could alter predic- 
tions of the critical temperature in an unknown way. 



Conclusion. Previous efforts have addressed the need 
of a curvature term in the liquid drop expansion of nu- 
clear masses, but no consistent interpretation was made. 
Some works state that it is unnecessary, and that it is 
enough to stop the expansion at the level of a surface 
term Q. Other studies give conflicting results, and even 
the sign of the curvature correction remains ambiguous 
0-01 • Some of these references do give results that agree 
with the ones here, but do not offer a physical picture. 

We demonstrate that the surface energy coefficient in 
the traditional liquid drop formula changes when differ- 
ent mass ranges are considered. The decreasing trend 
in the surface energy coefficient with increasing mass 
number is consistent with the presence of a curvature 
term. We present a consistent description of the curva- 
ture term's nature, determine its sign and demonstrate 
its presence in the nuclear masses. 

Simple physical arguments predict that the volume and 
surface energy coefficients should be equal. Without the 
introduction of the curvature term, the volume and sur- 
face energy coefficients appear to differ from each other. 
With the addition of the curvature term, the two coeffi- 
cients agree within error. 

The nature of the "Wigner" term linear with isospin is 
also considered. A slight change in the definition of the 
squared isospin, possibly quantum mechanical in nature, 
captures its relative magnitude without introducing an 
additional parameter. 

What is gained through these considerations is a 
streamlined physical picture of the liquid drop model. 
Consider the difference of the original liquid drop model 
in Eq. ([ij to the final equation presented here in Eq. (|11|) . 
Even though the latter appears more complicated, there 
are the same number of free fit parameters as the former. 
Instead of adding more and more terms to produce more 
and more exact representations of the nuclear masses, 
we have added a geometric physical picture and kept the 
same number of variables to obtain a more accurate re- 
sult. The lessons learned with this equation are more 
telling than letting all the parameters free. 
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